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In the last ten years several different types of actuators have been developed and fabricated 
for active control of rotors. This paper deals with a special hydraulic actuator system capable of 
generating high forces to rotating shafts via conventional bearings. The actively controlled hydraulic 
force actuator features an electrohydraulic servo valve which can produce amplitudes and forces at 
high frequencies necessary for influencing rotor vibrations. The mathematical description will be 
given in detail. The experimental results attained verify the theoretical model. Simulations already 
indicate the usefulness of this compact device for application to a real rotor system. 


1 Introduction 


Rotating machinery can be found in many fields of the industrial world; many could be im- 
proved with regard to speed of rotors/shafts, weight, noise, longevity and last , but not least, safety, 
by applying active vibration control. The key to success for active measures in the improvement ol 
the dynamics of rotating machinery lies in the availability of suitable actuators which have to satisfy 
the following requirements: the actuator must be capable of amplitudes in the range of the vibration 
amplitudes to be influenced on the one hand, and possess an appropriate frequency characteristic on 
the other hand. In the literature many papers deal with magnetic actuators, e.g. [3,4,5]. Magnetic 
actuators require a relatively large amount of space compared to the attainable forces. This deficiency 
may be avoided by the use of hydraulically controlled chambers. 

The main topic of this paper is the introduction of a newly developed hydraulic actuator which 
is able to apply forces to rotating shafts via the bearings. For producing forces (pressure variations) 
servo valves are used. The regulating distance is caused by elastic deformations of membranes (no 
cylinder piston arrangements). 

All influences which could be of practical relevance (oil compressibility, dynamic effects of 
the servo valve and the fluid itself) are taken into account in the mathematical description of the 
system. Applying the theory of plates, the strain in the membranes can be calculated. A limitation 
of the actuator system is given by the permitted strain which may not be exceeded in any area. 
Because of the actuator dynamics an appropriate enlargement of the equations of motion for the 
entire system is necessary. The equations of the actuator dynamics show that the chamber system 
already applies damping effects to the rotor system without any control input. With the aid of 
a well adapted controller the damping of the entire system can be essentially improved. This is 
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demonstrated by simulations. The input parameters, especially with regard to the actuator system, 
have been determined by experiment. 

For verifying the results achieved by computer simulations two different experimental test 
facilities have been constructed (see Fig. 4 and Fig. 5). One test rig contains the active chamber 
system. It allows the evaluation of its characteristic frequency. The experimental results which are 
attained verify the theoretical model. The second test facility consists of an elastic rotori structure 
and two different types of magnetic actuators. In further investigations the two test facilities will be 
combined in order to test the dynamic behavior of the elastic rotor controlled by the active chamber 
system. 


N omenclat ur e 


A 

system matrix 

A* 

membrane area 

b 

coefficient of the actuator input 

B 

control matrix 

c, d 

bearing coefficients 

C K 

lag coefficient of the control force 

cl 

membrane stiffness 

di 

passive damping coefficient of the actuator 

E 

identity matrix 

h 

force vector 

Epq 

leakage coefficient 

K v 

loop gain of the servo valve 

K* 

constant of the oil compressibility and the membrane buckling 

m,p,q 

system matrices of the second order system 

F* 

control force to the bearing housing 

Tl 

transfer function describing the fluid dynamics and the stiffness 

T P 

transfer function of the chamber pressure 

Tv 

transfer function considering the servo valve dynamics, the 


oil compressibility, and the membrane buckling 

T ualve 

transfer function describing the servo valve dynamics 

T x 

transfer function of the regulating distance 

Jl 

Jacobian matrix of translation of the bearing 

A P 

pressure drop in the valve 

AP S 

supply pressure difference 

A P v 

chamber pressure difference 

Q 

vector of generalized coordinates 

Q 

oil flow through the loaded servo valve 

Qv 

oil flow through the unloaded servo valve 

s 

Laplace-operator 

t 

time 

u 

control vector 
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U v 

control voltage to the servo valve 

U N 

nominal voltage of the servo valve 

z 

state vector 

x,y 

acting directions 

XLiVL 

regulating distance (bearing deflection) 

Uly 

natural frequency of the servo valve 

0 

rotor frequency 


servo valve damping coefficient 


2 Description of the Actuator Device 


The device of this compact system is shown in Fig. 1 by a schematic. It consists of two 
components: a commercially-available servo valve, which is able to transform electrical signals into 
flow- and pressure variations of the hydraulic fluid, and an elastic chamber system, which transduces 
the pressure variations into displacements and into forces acting on the outer bearing housing (in 
Fig. 1 only one control direction is displayed). 

The servo valve is supplied by the system pressure which is to be kept constant, and the return 
pressure R & 0. If the spool of the valve is in .themeutr all position' (Fig. 1} the pressure in the chambers is 
half of the system pressure. After a shift of the spool by an input voltage Uv ^ O.a pressure difference 
appears between the two output orifices. Because of the elasticity of the membranes this enables a 
shifting of the rotor in the radial direction. The chamber system itself consists of four cylindrical 
chambers which are equally spaced in a circle around the bearing housing. Each chamber is sealed on 
the top and bottom by an elastic membrane. In order to decouple forces into both acting directions, 
the bearing housing is supported against the membrane system by linear roller guides. The influence 
of friction is thereby reduced as well. Alternative support solutions are given by design variations. To 
eliminate the linear roller guides the bearing can be mounted via elastic rods to the outer housing. 
The two opposing chambers are both controlled by one servo valve (one valve per force direction). 

This compact actuator system is capable of generating very large forces and can thereby influ- 
ence even large turbines weighing several tons. Actuator pistons as an alternative to the deformable 
chambers have the drawback that they possess relatively large moving masses. In addition, sealing 
problems and friction forces on the contact surfaces would occur. All these influences have negative 
impacts on the frequency characteristic. An appropriate design of the membranes enables the desired 
amplitudes of motion without exceeding the stress limits of the material. 


3 Theory 

With the aid of the theory of plates the calculation of the stress in the membranes caused by 
a given load can be performed, e.g. [8]. Experiments have shown that the real stresses are always 
lower than the calculated stresses. This is caused by the elastic fixing of the membrane in the real 
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system compared to the theoretically stiff fixing. For calculating the stresses the load must be known. 
Therefore it is necessary to describe the dynamics of this active machine element. 

3.1 Chamber System Dynamics 

3.1.1 Servo Valve 


The dynamics of the servo valves are given by the manufacturer (MOOG, see [2]) as a PTV 
element. The valve flow Qv for the unloaded valve can be given as a function of the input voltage 
Uy in the Laplace- domain (with s as a complex variable) 


Qv — ffuT i valve ' Pv ) *P' valve 


1 + 2k s + 


( 1 ) 


where I\ v is the loop gain of the valve, which depends on the supply pressure A Ps = P — R (see 
Fig. 1) and the geometry of the valve itself, lo v is the natural frequency of the valve, and £„ is the 
damping ratio of the valve. In the case where the valve is loaded, the fluid flow depends on the output 
pressure A Py (load differential pressure). For this reason the equation for the fluid flow passing a 
pressure- regulating valve (here between the spool of the valve, inlet- and outlet orifice, see Fig. 1) is 
applied. In general the fluid flow Q is proportional to the square root of the differential pressure A P 
[1], 


Q ~ VAP 


(2) 


Applying the relation that the sum of the pressure drops in the valve A P and the load pressure 
A Py is equal to the drop of the supply pressure A Ps, 

AP+AP V = AP S , ’ (3) 

and considering eq. 2 and the fluid flow properties in the valve we obtain 

Q ° Qv '\f^KFs ' (4) 

After a linearization in the neighborhood of A Py = 0,eq. 4 becomes 

Q = Q v — K pq • A Py , K pq leakage coefficient . (5) 

The linearization applied in eq. 5 is valid for A Py < 5A.P5 w ^ ere K vq theoretically a 
function of Q v . For low input voltage Uv < JqUn {Un = nominal voltage) this dependance can be 
neglected. Fig. 2 shows the course of the fluid flow Q as a function of Qv and the differential pressure 
of the valve A Py as given by the manufacturer. The slope of these curves at the operating point 
(A Py — 0) is given by K vq appearing in eq. 5. In contrast to the nonlinear eq. 4 in practice one can 
establish that even at Uv = 0 (which results in Qv = 0, eq. 1 ) the fluid flow is non-zero for APv 7 ^ 0. 
This is caused by leakage where the spool is in the neutral position. The linear eq. 5 takes this into 
account. Fig. 2 shows this behavior of the valve by the line through the origin with a slope which is 
non-zero. For derivation of the equation of motion and the transfer characteristic of the entire system 
(see later), the linear relation by eq. 5 is used. The nonlinearity of the relation between the fluid flow 
Q and the pressure drop A Py can be taken into account by K pq as a function of Qv- 
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3.1.2 Influence of Fluid Flow, Membrane Stiffness and Oil Compressibility 

The application of the continuity equation results in an additional relation for the fluid flow 

Q, 

Q = A* ■ x L + I<* ■ APv . (6) 

A* is the characteristic membrane area, and K* represents the oil compressibility and the deforma- 
tion of the membrane into a direction that doesn’t lead to a movement of the bearing (membrane 
buckling), see Fig. 3. K* is dependent on the oil volume Vl, the bulk modulus of the fluid 0k, and 
the membrane geometry characterized by Km, 

K * = 77T + K m ■ ( 7 ) 

4 0k 

Simulations have shown that the fluid losses caused by oil inertia, friction and losses at the inlet 
and outlet orifices are negligible. So the pressure difference at the outlet orifices of the valve is 
approximately equal to the chamber pressure difference, constant A Py. 


,3.1.3 Recapitulation of Equations 

The equilibrium of forces applied with respect to the membranes (entire arrangement of two 
membranes acting in one direction) leads to 

F* = A* ■ AP V -c L -x L . , (S) 


where cl is the entire membrane stiffnessintheforce'directioniand F 1 * the control force to the bearing 
housing, see Fig. 3. Including this equation we obtain F* as a function of the control voltage Uv and 
the membrane deflection, respectively; the bearing displacement xl in the; direction of the force 

F* = Fv-Uv - T L -x L . (9) 


Tv is the transfer function describing the dependence on the control voltage 

Ty = T valve 


A* K v 
Kva 


1 + Ck • s 


( 10 ) 


and Tl is the transfer function considering 
the arrangement 


T l = 


the influences of the fluid dynamics and the stiffness of 

+ C L - ( 11 ) 


Kr 


1 + Ck ■ s 


The constant Ck is l'epresented by the 
membrane buckling K* eq. 7, 


leakage coefficient 


Ck — 


K* 

K Pq 


K pq 


and by the oil compressibility and 


( 12 ) 
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3.1.4 Representation in Time Domain 


For including the actuator system into the complete system, the representation of eq. 9 in 
the time-domain is more convenient. With the equations above we get a set of differential equations 
which describe the control force F* as a function of the input voltage Uv and bearing displacement 
XL, 


F* = A*AP v -c l x l , 

ICAPv + K pq AP v + A*i L = Q v ( 13 ) 

— ^ Qv + Qv + Qv = K v llv 

u v 

The coefficient cl describes the additional stiffness effects with regard to the entire system (see later). 

3.1.5 Simplification and Interpretation of Equations 

For designing control concepts and to explain the fundamental operation, the equations will 
be simplified. In order to obtain a simple relationship between the control voltage Uv, the regulating 
distance xl, and the control force F* the PT 2 -behavior of the servo valve will be approximated by a 
linear relationship (F va i ve = 1)- This is valid at frequencies under the cut-off frequency of the valve. 
This simplification leads to the following form 

CkF* + F* = b • Uv — cl ■ xl — <1l‘Xl ( 14 ) 

where b = -i, Kv and d,L = 4— + CkCl ■ 

J\pq JApq 

Tb' 3 constant Ck representing the oil compressibility and the membrane buckling in the equa- 
tion above has an important influence. If Ck is negligible the actuator works as a passive spring and 
damper unit characterized by the coefficients cl and di- Additionally a control force Fc acts that 
is proportional to the input voltage Uv (see Fig. 9). In the other case, if the frequency uk = ^ 
is below the frequency range of application, the actuator possesses an integral behavior. As a result 
the feedback of velocities and accelerations leads to stiffness and damping effects. In this case or 
if t ok is inside the range of application (PTi behavior), a proportional behavior (see Fig. 9) can be 
achieved by an internal phase shifting controller (e.g. feedback of pressure APy ) or an additional 
passive throttle between the two output orifices. 

By an appropriate design of the membrane (constants cl , A *, Km), of the oil volume Vl, of 
the supply pressure APs , and by the choice of the suitable valve size (constant I( v and Kpq), the 
actuator can be tuned in an optimal way. In the case of a sudden pressure drop ( APs — » 0) the rotor 
bearing is always sufficiently supported by an adequate membrane stiffness cl- 


3.1.6 Transfer Characteristic of Actuator 

In the following the transfer characteristic of the complete actuator system is investigated 
applying a sinusoidal input voltage Uv to the valve. Thei purpose ofithisisto get a comparison to 


274 



the experiments which are described in chapter 3.1.6 and 5. The resulting outputs are the regulating 
distance xl and the pressure difference A Py between the opposing chambers. Introducing a mass 
which simulates the mass of the rotor, including the support system, the control force appearing in 
eq. 9 has to be equal to 

F* = m xl • s 2 . (15) 


The transfer functions for the regulating distance F x and the chamber pressure Fp can be expressed 

b • F i naive 


j- _ XL _ 

Uv Cfcm • s 3 + m ■ s 2 + d,L • s + cl 


(16) 


F P = 


APy m ■ s 2 + cl 


U v A * 

The coefficients b and d,L are given in eq. 14. 


F x 


(17) 


4 Test Facilities 

Fig. 4 shows the test rig which was designed to investigate the chamber system experimentally. 
The investigation is conducted only in the vertical direction. The facility mainly consists of the 
housing 1 where the four oil chambers 2 are integrated. The distance rod 3 performs the function of 
the bearing housing and assures that the opposing chambers both move the same absolute distance. 
For investigating the influence of masses additional masses can be mounted on the distance rod 3. 
The influence of friction caused by the linear roller guides (see Fig. 1) can be studied by loading 
the chambers for the horizontal direction with an appropriate pressure in order to press the linear 
roller guides via distance sleeves 4 on the distance rod. The entire arrangement is mounted on 
the foundation including the servo valve 5. Servo valve and chambers are connected via pipes 6. 
Between actuator and hydraulic control unit (not shown in Fig. 4) are the connections to the first 
stage pressure unit, main control pressure unit (pressure P shown in Fig. 1), return pressure unit 
(indicated by pressure R in Fig. 1), and an additional pressure supply for the linear guiding bearings. 

To enable a broad range of experimental investigations, variations of the supply pressure A P5, 
membrane thickness h, contact pressure for the guiding bearings, additional masses m (see eq. 15) 
and of course the control voltage Uv (function generator) can be conducted. 

The measurement equipment installed allows one to gather the following information: regula- 
ting distance with the aid of an inductive displacement sensor 7, and oil pressure in the chambers 
by piezo electric pressure sensors 8. Additional information is given by acceleration pick ups (not 
shown) mounted on the distance rod and by strain gauges to determine the strain of the membranes. 

A special hydraulic unit was designed to assure that the supply pressure is constant. This is 
required because the fluid flow changes with frequencies up to 200 Hz and higher. These requirements 
can be met with the aid of a controllable radial piston pump and an: appropriate reservoir. 

The test rig shown in Fig. 5 is used to investigate different control concepts. It consists of an 
elastic rotor structure 9, a magnetic bearing 10 for simulating different excitation forces, the bearing 
unit 11, and actuator system which is presently realized by electromagnetic actuators 12. A further 
description of the test rig and the test results attained by this facility are! given e.g. in [3]. 
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The active chamber system will be integrated in the test facility Fig. 5 after conclusion of the 
experiments with the test rig Fig. 4. For this purpose the chambers for the horizontal direction have 
to be activated by a second servo valve, and the distance rod has to be replaced by the bearing unit. 
Finally the electromagnetic actuators 12 are substituted by the active chamber system. 


5 Transfer Characteristic 

In the following the transfer functions J~ x and J-p.e q. 16 and 17, for the regulating distance xp 
and the chamber pressure A Py will be investigated. The goal is to compare theoretical and experi- 
mental behavior of the actuator and to determine the system parameters. These are the membrane 
stiffness Cp, the leakage coefficient K pq , the loop gain of the servo valve K v and the coefficient I\* 
which describes the oil compressibility and the membrane buckling. The simulation is based on the 
data shown in Fig. 6 which are equivalent to the data of the test facility. The results can be seen in 
Fig. 7. 


5.1 Theoretical Results 


In the plot’ of the amplitude of the regulating distance in Fig. 7, one can recognize a PTj 
behavior up to a frequency of about 170 Hz. Its cut-off frequency is very low (about 2.2 Hz) and is 
given by see eq. 14. Beyond this frequency an extremum occurs in the amplitude at approximately 
300 Hz. It corresponds to a pole in the transfer function of xp and APy. This extremum results from 
a function a(s ) which can be obtained by the denominator of eq. 16. This function has the value 
e*(s) = 0 at the frequency of u> a = 2n • 287 Hz (s = iu>), 


, . o m dp 

c*(s) = m-s + 77- • * + — 
wr Wv 




0 


dp 1 

CWm “ 4 C\ 



(18) 

(19) 


Eq. 19 shows that this frequency is mainly dependent on the membrane area A*, on the mass m 
(simulating the rotor mass including the bearing unit) and on the oil compressibility and membrane 
buckling K*. The servo valve has no influence. The expression tU- can be interpreted as a spring 
coefficient of the oil support and theunembranebuckling. This natural frequency causes a phase shift 
of 180°. The rest of the. curve is characterized by a drop in the amplitude due to the PT 2 -behavior 
of the servo valve, see chapter 3.1.6. In the plot of the chamber pressure a further zero point 
can be seen at 110 Hz. It corresponds to the natural frequency of the one-DOF-oscillator, function 
(3{s) from eq. 17 

/3(s) = m ■ s 2 + cp = 0 (20) 

( 21 ) 

This characteristic frequency has no effect on the regulating distance and on the controllability of 
the rotor system. 
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5.2 Comparison of Theory and Experiment 

Fig. 8 shows the measurements of the regulating distance xl and chamber pressure A Py. The 
results are obtained by noise excitation for the input voltage Uv- 

A comparison with the results attained by theory shows good agreement. The characteristic 
frequencies at 110 Hz and 293 Hz are due to the effects described above (functions a(s) and (3(s), 
eq. 19 and 21). This shows that all important influences on the transfer characteristic are included 
in the transfer functions eq. 16 and 17. The step shaped curve at lower frequencies is only caused 
by the discrete evaluation of the output signals by the measurement electronics. 

Measurements of the transfer characteristic of the servo valve (described by a PT 2 -system) 
showed a better amplitude behavior at frequencies beyond the cut-off frequency than predicted by 
the PT 2 -behavior. 

The experiments have shown that the simplified equations for describing the performance of 
the active chamber system, eq. 14, are valid up to the cut-off frequency of the servo valve (here 
about 200 Hz), see chapter 3.1.5. For a higher frequency range it is necessary to use a servo valve 
with a higher cut-off frequency and to minimize the oil compressibility and the membrane buckling 
by minimizing the oil volume and by an appropriate design of the membrane. 

6 Example of Application of the Actuator 

6.1 Formulation of the Complete System 

In designing a controller, a complete mathematical formulation of the system is 
assumed. A very efficient method of describing the rotor system; to be controlled is modeling it as a 
hybrid multibody system (see [6]). This leads to the second order differential equation 

Mq(t) + Pq(t) + Qq(t) = Y^hi(q,q,q,Cl,t) (22) 

i 

with the vector q e IR/ of the generalized minimal coordinates (/ is the number of degrees of freedom), 
M e IR.™ is the mass matrix, P e IR/’^ is the matrix of velocity-proportional forces, Q e IR/’^ is the 
matrix taking into account the displacement-proportional forces, hi e DR/ is the vector of the i-th 
external force (e.g. control force), Cl is the rotor frequency and t indicates the time. 

The control forces acting indirectly (via the rigid bearing unit) on the rotor can be expressed 

by 

h B = J T L - F i , (23) 

. r y . 

where J l £ IR 2 ’^ is the Jacobian matrix of translation belonging to the subsystem bearing unit and 
F* and F* the control forces given by differential eq. l3.acting in the x- and y-directions . The radial 
movement of the bearing described by x B and yi (perpendicular to x B ) can be formulated with the 


111 



aid of the Jacobian-matrix as 


XL 

VL 


Jl • q 


(24) 


Introducing an enlarged state space vector 


z = 


q j Qvxi Qvyi q lQvxiQvyi^Pvxi^Pvy 


z e 


IR 2 /+6 


(25) 


with Qvx,Qvy&s the valve flow serving the x- and', y-directionis , respectively (eq. l),and A Py x , APv y as 
the valve output pressure supplying the x- and y-directions, respectively , weobtain the state equation 


z = Az + Bu 


(26) 


where A represents the system matrix, B the control matrix and 

u = [U Vx ,U Vy ] T (27) 

the control vector consisting of the control voltages Uv x proper for x- and Uv y proper for the in- 
direction. 


The matrices used in eq. 26 can be expressed as follows 



0 

0 

Ef 

0 

0 


0 

0 

0 

E 2 

0 

A = 

-M~ l Q 

0 

-M~ l P 

0 

M~ l J T L A 


0 

-u 2 v e 2 

0 

— 2£ v u> v E2 

0 


0 


A* r 
K* J L 

0 

-§*^2 


0 

0 


K v u 2 v E 2 

0 

where the submatrix Q is 

Q = Q + J Tj L ’ C L 


(28) 


(29) 


(30) 


cl is the membrane stiffness appearing in eq. 13 and M, P, Q are the system matrices from eq. 22. 
The matrix Ej e IR JJ represents the identity matrix. 


6.2 Amplitude Characteristic of the Rotor 


A simple rotor running in journal bearings is used. The rotor is supported by the active 
chamber system on its left side (Fig. 9). In the modeling of the rotor system. mass-, stiffness-, and 
gyroscopic effects and additional damping by the journal bearings are considered. As admissible 
shape functions only two mode shapes are taken into account. In most cases higher modes of real 
rotor systems are damped very well because of bearings or material damping. 
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In Fig. 10 are plotted frequency response functions versus rotor frequency. The exciting force 
is caused by an unbalance of the rotor. In Fig. 10a the scaled amplitudes ■— at the runner position are 
displayed without active forces (solid line) and with active forces (dotted line). Fig. 10b 
shows the scaled nutation angles — (slope of the elastic line) at the same location. The displacement 
amplitude shows a significant peak at the first natural frequency (bending vibration) and at the 
second natural frequency (nutation, indicated by the nutation angle). The feedback is realized by 
displacement, velocity, and pressure signals taken only at the actuator location. The controller was 
designed as a ’’constant controller” for the rotor frequency 100 Hz. Using this very simple controller, 
the results already (demonstrate that a considerable reduction of the resonance amplitudes can be 
achieved. 


7 Conclusion Remarks 

The key to attaining the desired influence on the dynamics of rotor systems lies in the availa- 
bility of suitable actuators. The development of the active chamber systems introduced in this paper 
will be an essential step towards industrial applications. This actuator works as a spring (membrane 
stiffness) and damper (hydraulic effects in the servo valve), and in addition to these passive, forces a 
control force is acting proportional to the input voltage at the servo valve. One of the main goals of 
the paper is to confirm the theoretical model of this new actuator concept by experiment. The theore- 
tical and experimental results achieved by these investigations show a very good conformity. Further 
experimental investigations of the actuator system included in the rotor system will demonstrate the 
usefulness of applying this active machine element for active vibration control of rotors. The use of 
a more sophisticated controller (additional measurement information, adaptive controller realized by 
special electronic equipment, hybrid controller) will show a wide range of possible applications for 
such an actuator device. 
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Figure 1: Active chamber system with servo valve 
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Figure 5: Test rig for rotor control 
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Figure 6: Data for the transfer characteristic of the actuator 
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Figure 9: Geometry and mechanical model of the investigated active rotor system 
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Figure 10: Response of the unbalanced rotor running through the critical speeds 
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